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1 Introduction 


Today it is clear for many physicists that one of the main goals of the future CERN Large 
Hadron Collider {LHC) is to hnd as much information as possible about the nature of 
the Standard Model {SM) Electroweak Symmetry Breaking {ESB) [|l|. In spite of the 
huge amount of data obtained in the last years at the Large Electron-Positron Collider 
{LEP) it is very few what we really know about the ESB. The proposed mechanisms 
include ideas such as supersymmetry (see and references therein), technicolor and 
many others. Therefore it would be very interesting to have some model independent 
framework to make a phenomenological description of the ESB physics. 

In fact such a framework exists, at least for the strongly interacting case i.e., when no 
light modes are present in the ESB. It has been developed in the last years and used 
to describe the scattering of the longitudinal components {LC) of the electroweak gauge 
bosons m as well as the precision tests of the Standard Model coming from LEP |^. It is 
based on the application of the chiral Lagrangians or Chiral Perturbation Theory {xPT 
) 0 (previously invented for the description of the low-energy hadron interactions) to the 
dynamics of the ESB Goldstone bosons {GB). 

One assumes that there is a physical system with a global symmetry group G which is 
spontaneously broken to some other group H. This global symmetry breaking drives the 
gauge electroweak SU{2)l x U{1)y symmetry breaking down to U{l)em through the well 
known Higgs mechanism. The only election for the G and H groups compatible with the 
presence of the SU{2)i^ji custodial symmetry (to have a p parameter naturally close 
to one) and the existence of three massive gauge bosons after the ESB i.e., the , W~ 
and ZO, is G = SU{2)l x SU{2)r, H = SU{2)l+r. 

Then the low energy dynamics of the ESB is described by a SU{2)l x U{1)y gauged 
non-linear sigma model {GNLSM) including an arbitrary large number of terms in the 
action with different number of derivatives of the GB helds and electroweak gauge bosons. 
The corresponding couplings (parameters) encode the dynamics of the ESB sector of the 
SM and must be renormalized to absorb divergences. However, at low enough energies 
only a small number of terms (and couplings) are needed to eliminate all the divergences. 
In principle the values of these couplings or parameters could be obtained from the un¬ 
derlying theory or directly by htting them from future experiments. 

In this work we will study the possibilities of the LHC for measuring these param¬ 
eters and the expected corresponding errors. In particular we will concentrate in events 
producing Z°Z° or pairs. As it was mentioned above, the application of the chiral 

Lagrangian technique to the production of electroweak gauge bosons is not new and has 
already been considered at the literature 0. However, all the applications worked out 
until now are based in the so called Equivalence Theorem {ET) ^j. This theorem relates 
the S matrix elements of processes containing electroweak gauge bosons LC with the 
corresponding processes with GB. However, in a recent work concerning the formu¬ 
lation of the ET in the context of xPT-i the Equivalence Theorem is severely restricted to 
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a narrow energy applicability window. It can also be applied in the high energy domain 
together with the yPT but using some non-perturbative technique, like dispersion rela¬ 
tions 1^ or the large N limit []^, in order to have an appropriate unitarity behaviour of 


the amplitudes. 

For this reason we consider prioritary to apply directly the chiral Lagrangian descrip¬ 
tion of the ESB without using the ET. The main problems of this approach are two: 
First one has to include explicitly the gauge degrees of freedom in the model which makes 
the computations extremely more difficult. Second, one has to restrict the results to the 
low energy region where standard xPT can safely be applied thus losing many higher 
energy events. The advantage is that the values of the htted parameters will be more 
reliable since one is not using the ET complemented with some non-perturbative method. 

The plan of this work goes as follows: In section 2 we introduce the chiral effective 
lagrangian to be used with its parameters and the corresponding Feynman rules. In section 
3 we compute the cross-section of the LHC subprocesses that are relevant for measuring 
the parameters. In section 4 we consider the signatures and the possible backgrounds. In 
section 5 we show how we compute the total number of expected events at the LHC from 
the subprocess cross-sections for the signals and backgrounds. In section 6 we discuss the 
sensitivity of the machine to the different parameters, we dehne the optimal cuts and we 
compute the statistical signihcance of the different parameter measures and estimate the 
statistical errors. In section 7 we consider other sources of systematical errors such as the 
uncertainty on the proton structure functions and in section 8 the effect of the running 
of the parameters. Finally, in section 9 we review the main conclussions of our work. 


2 The effective Lagrangian and the Feynman rules 

In this section we consider the most general effective Lagrangian for the ESB compatible 
with the SU{2)l x U{1)y local symmetry of the SM and the breaking pattern SU{2)l x 
U{1)y —^ f^(l)em- The well known chiral Lagrangian describing the ESB of the SM, can 
be written as an inhnite expansion with terms of increasing number of gauge helds and 
derivatives of the Goldstone bosons {CB), with an inhnite number of arbitrary parameters. 
This chiral Lagrangian can be seen as a low momentum expansion for the corresponding 
Green functions. At some given order in the number of GB derivatives one can work 
only with a hnite number of terms and parameters. In this case the model can only 
be applied to much smaller energies than Attv which is the parameter controlling this 
expansion (dvru ~ 3TeV since v ~ 250GeV). 

The hrst term of this effective Lagrangian (0(p^)) is given by: 

£( 2 ) = ^tr[D,U{D>^Uy] ( 1 ) 
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and by the Yang-Mills Lagrangian: 


C.YM = - \Tr\F^,F^''\ 


( 2 ) 


We choose the GB parametrization as an unitary matrix IJ belonging to the quotient 
space ^17(2)i x S\J{F)nl 


T .71 

U = exp(i-) 

V 

where r = (Ti,r 2 ,r 3 ) are the Pauli matrices and tt = (7ri(a;),7r2(a;),7r3(a;)) represents 
the triplet of GB. 

The covariant derivative D^U is dehned as: 


D^U = d^U + tgW^U-tg'UY^ 

Here, and are the SU(2)l and the U{1)y gauge fields given by: 


W^ = 


W^f 

2 


Y — 

^ 2 

where IT^ = (hP^, W‘^, W^) represents the triplet of SU{2)l gauge helds. 
As usual, the covariant held strength tensors are dehned as 


F^,{x) = d^W,{x)-d,W^{x)+ig[W^{x),W,{x)] 

Bi,y{x) = di,By{x) - dyBf,{x) 

The transformation properties under the SU{2)l x U{1)y gauge transformations are 

= gL{x)igWf,gl{x) + gL{x)d^gl{x) 
ig'yl, = ig'y^l + 9Y{x)d^g\r{x) 

U' = gL{x)UgY\x) 

where 
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gdx) 

9y{x) 


_ gi6»fe(a;)r*/2 

_ gi/3(x).T3/2 


Note that, as expected, the GB helds, tt*, transform non-linearly . 

In order to construct the chiral lagrangian to order O(p^) we dehne, following Longhi- 
tano in [T^, the quantities T, and V^O{x)\ 


T = Ut^U^ 

= D^UU^ 

V^O{x) = d^O{x)+ig[W^{x),0{x)] 

Then the complete electroweak chiral Lagrangian with the whole set of SU{2) x U{1)y, 
Lorentz, (7, P and T invariant operators up to dimension four has the form: 

C = + Cym + 1 + X/ (^) 

i=l 

In this equation, the different terms Ci are the SU{2)l x U{1)y invariant functions of 
the gauge vector bosons and the GB fields, containing four derivatives, whereas C[ has 
dimension two. They have the following expressions: 

C\ = ^g‘^aov^[Tr{TVf,)f 

Cl = ^-g‘^aiB^,Tr{TF^’^) 

C 2 = Uga2B^,Tr{T[V^^,V'^]) 

£3 = iga:iTr{F,,[V^^,V'']) 

A = «4[Tr(AK)]' 

£5 = a^[Tr{y^V^)f 

A = a6Tr[(AK)]Tr(TV^^)Tr(Tf"") 

Cj = ajTr[{V^V^)][Tr{TV'^)f 

A = \g^a,[Tr{TF^,)f 
A = ^*^/a9Tr(TA.)Tr(T[f"^n) 
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Cio = ^«io[Tr(TV;)Tr(TK)]' 

£11 = a^iTr[{V^V^f] 

C12 = ^a, 2 Tr{TV^V,V^)Tr{TVn 

Cis = ^a,3[Tr{TV^Vn? (4) 


As it is well known there is some arbitrariness in the choice of a particular base of 
invariants. Using the equations of motion it would be possible to eliminate the operators 
Cii, C 12 and Ci 3 redehning the rest of the terms (see for example [|I4|). 

So, a set of ak {k = 0, ..13) parameters have appeared in the definition of this effective 
Lagrangian. We display in Table |l| the relations between some different sets of the chiral 
Lagrangian parameters used in the literature, i.e. by Longhitano [|T^], Feruglio |T^ and 
Dobado et.al. p. 

The theory defined by the Lagrangian above is non-renormalizable in the strict sense, 
but the divergences appearing when we calculate at one-loop with the chiral Lagrangian 
at lowest order eq. (0) have the same form that some of the terms obtained at tree level 
in the following order in the chiral expansion. The counterterms needed to reabsorb the 
divergences generated to the one-loop level with in the Landau gauge were obtained 
by the authors in JT^. They have the same form that C[, Ci, £, 2 , £ 3 , £a and £ 5 , whereas 
the rest of the terms in eq.(^ are not needed for renormalization. 

We will describe the interactions between gauge bosons and fermions by the same 
Lagrangian than that of the S'M, which can also describe the couplings between fermions 
and scalars by means of the usual Yukawa terms. However the latter are not used in our 
calculations since we work all the time with massless fermions. 

Once we have set the Lagrangian to be used, the next step is to define the quantum 
theory. This can be done in a standard way using the Faddeev-Popov method and choosing 
some appropriate covariant gauge. The gauge fixing functions are constructed by 

requiring that the terms of the form HT 9 ^ 7 rj(x) and B^d^'K 3 cancel out; thus yielding: 


z = l,2,3 

Now we can write the gauge-fixing and Faddeev-Popov terms as: 


Cgf = 

i=i 


( 5 ) 
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m, n = 0,1, 2, 3 


^FPG J d yC m\X)-^ 

Therefore, the effective action and the effective Lagrangian (non-linear SM Lagrangian, 
d^NLSu) take the form: 

Seff = J d^a;[£(a;) -|- Cgf + ^Ifpg] (6) 

C-NLSM = Cp‘'’ + l + + d^YM + d^GF + C-FPG + (7) 

In the following we will to work in the Landau gauge (^vv = = 0)- la this gauge the 

perturbative tt propagator remains massless and the Faddeev-Popov ghosts only couple 
to the gauge held. There is not direct coupling of the tt held to the ghosts. Moreover, the 
counterterms necessary to cancel divergences with at one loop are gauge invariant 
functions of the gauge helds and the GB alone (the £i, £, 2 , £ 3 , £ 4 , and £5 terms 
mentioned above). In other gauges, or other parametrizations of the coset space, the 
counterterms could also be functions of the ghost helds, and their structure would be 
determined using the Becchi-Rouet-Stora {BRS) invariance. 

Thus, the gauge hxing and the Faddeev-Popov term in the Landau gauge read: 

£fpg = - geijkB^c’d^W^'^ + df,cld^c° (8) 

In the Appendix we show the Feynman rules derived from the Lagrangian £nlsm 
eq.(0) that we use for our calculations. They correspond to Figures SI and As usual 
it has been useful to redehne the gauge helds as: 


Z, 

A, 


cosOwW^ — sinOwB^ 
+ cosOy/B^ 


where the weak angle is dehned by 


tgBw 


I 

9 
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and the GB as: 


TT 


± 


TT 


0 


TT^ =F 


3 The subprocesses cross-sections 

In this work we are interested in studying of the different cross sections that contribute 
to the final states and W^Z^ at pp colliders such as LHC since they are the most 

promising from the experimental point of view. In both cases we focus our attention on 
the so-called gold-plated events, where the produced bosons decay to the leptonic hnal 
states I = e,fi, v. In spite of the very small branching ratios {BR), 


BR{Z^Z^ lll'l') = 0.0044 

BRiW^Z^ = 0.013 

these are the most interesting events since they are much easier to detect than the 
hadronic channels. 

On the other hand, we will assume that it is not possible to detect experimentally the 
different polarizations of the gauge bosons in the hnal state. Therefore we have to include 
in our computations all the Z^ or gauge bosons, either transversally or longitudinally 
polarized. As it was mentioned in the introduction, instead of using the ET to calculate 
the collisions of the weak bosons, as it is usually done, we have computed the amplitudes 
with gauge bosons as initial and hnal states and then projected them in all their polar¬ 
izations. Thus, as we are not going to apply the ET we only consider the maximal bound 
on the energy required for the applicability of the chiral Lagrangian formalism, but we do 
not have any low energy bound. In this way we can take into account the inhuence of the 
chiral lagrangian parameters in all the gauge bosons polarization channels and not only 
in the longitudinal ones as it is usually the case. The main disadvantage in our procedure 
is that the calculations become much more complicated. This fact leads us to work only 
at tree level to make the computation more accessible. 

In principle, according to the spirit of yPT, the one-loop corrections coming from 
the lowest order lagrangian should also be included since they are O(p^). However, the 
computation of these corrections is extremely involved. In addition, they have not de¬ 
pendence on the chiral parameters. As it will become clear later, our main interest is to 
study the dependence of the number of events on the chiral parameters in order to see 
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which of them could be measured at the LHC. For this reason we do not expect that our 
results concerning the measurable parameters would change too much when the one-loop 
corrections are included but, of course, the precise number of events will do. 

Later we will also study the effect of the running of the parameters on our results 
which in some way takes into account part of the one-loop contribution, but, in any case, 
our approach is the simplest one including the right dependence on the chiral lagrangian 
parameters. 

Now we describe how we have calculated the cross sections of the above mentioned 
processes. From the quantized chiral Lagrangian that appears in eq.(|^ we obtain the 
corresponding Feynman rules (see Appendix). Then we compute, at tree level, the ampli¬ 
tudes for the different subprocesses contributing to the or W^Z^ hnal state, using 

a REDUCE code, and then the corresponding cross-sections. We calculate them in the 
center of mass frame, with a total energy -\/I. There, the four momenta of the interacting 
particles are: 


Vi = (^i,0,0,p) 

P2 = {E2,0,0,p) 

P 3 = {E3,psm9,0,pcos9) 

P 4 = (E 4 , —p'sin^, 0, —p'cos9) 


Here, Ei are the particle energies, and p and p' are, respectively, the magnitudes of 
the three-momenta of the initial and the hnal states. 


p 


7712^) + {m2^ - 



p 


- 2s(m3^ m^) + {rrii - 


S = (p^+p2)^ 

i = (pi - p^Y 

U = (pi - P4)2 


where s, t and u are the Mandelstam variables. The differential cross section is given 
by this expression: 



( 9 ) 
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where M is the helicity amplitude, and the S symbol refers to the sum of all the dual 
polarizations and the average on the initial ones. In the case of electroweak gauge bosons 
scattering, where we separate the transversal and longitudinal polarizations in the initial 
state (as we will see in section 5), we take into account the different helicities contributing 
to TT, TL, LT and LL polarizations of the initial bosons, to do this average. We also 
include the corresponding branching ratios to the gold-plated events. 

Let us consider now the processes corresponding to the elastic collision of gauge bosons 
(fusion processes). Symbolically this reaction can be written as: 

-^1)^2 (P2; -^2) \‘ i ) V /^{ pi , A4) 

We choose the following polarization vectors in the helicity base (Aj = ±1, 0): 

4 = ^(0,Tl,-i,0) 

e? = ^(P, 0 , 0 , El) 

mvi 

4 = ^( 0 ,Tl,i, 0 ) 

4 = -^(p,o,o,-E2) 
mv2 

€3 = ^^( 0 , =fcos 6 ', —i, ±sin 6 *) 

v 2 

Co = -(p', E3sin6*, 0, E3COS6*) 

mvg 

I 

64 = ^^(0, =fcos 6*, i, ±sin6*) 

V 2 

X 

64 = - (p', —E4sin6*, 0, —E4COS6') 

mv4 

The number of independent helicity amplitudes in fusion processes of massive gauge 
bosons, such as Z°Z°, W+W- Z°Z°, W+Z^ W+Z^ or IT-ZO ^ W-Z^, 

is 3^ = 81. The photon has only two polarized transverse states and therefore, in the 
processes W^Z^ or W~'y W~Z^ the number of helicity amplitudes is 2 x 

3^ = 54. As our chiral Lagrangian affecting the boson scattering is invariant under 
(7, P and T transformations, we can derive many relations between different helicity 
amplitudes. However we have not used those relations to save computations. Instead we 
have calculated all the helicity amplitudes in each process and the relations between them 
have been only used to check our results. 
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3.1 final state 

Now we concentrate in the different processes that contribute to the final state 


99 

qq 


7^0 ^0 
Zj 'Y' ^ 'Y 

ryO ryO 

Zjt^L 

T'O 70 

zlzl 

W^WE 

WtWE 

W+WE 




Note that we separate the contribution of the different polarizations channels in the 
initial state, since different polarizations will have different luminosities in pp collisions. 


In particular we use the Weizsaker-Williams and the effective W approximation 
to compute the Z 1 V 2 pair luminosity in the LHC beams. Therefore, we have to divide 
these cross-sections into the contributions coming from the different gauge boson initial 
polarizations. 

On the other hand, it is well known that at supercollider energies the one-loop process 
gg ^ ZZ is not negligible. The gluon-gluon fusion cross section was calculated in the 
Minimal Standard Model, MSM (with just one Higgs doublet) via one-loop of quarks 
by Glover and Van deer Bij 0. The corresponding diagrams are shown in Figure 


As we are using a chiral Lagrangian description of the ESB sector we do not include the 
contribution coming from diagrams with the Higgs boson in our computations. 

The quark-antiquark annihilation represents the main source of Z^Z^ pairs in pp col¬ 
liders like the LHC. As in the case of gluon fusion, there is no dependence on the chiral 
parameters in the calculation, at our level of approximation. The cross section at tree 
level, which only receives contribution from the t and u channels (Figure has the same 
well-known expression as in the MSM. 

The Z^Z^ fusion calculated using the chiral Lagrangian, only receives contribution 
from the vertex Z^Z^Z^Z^ (Figures ^ and |^). The only dependence on the effective 
Lagrangian parameters is through the following combination eq.( 


0:4 -f- tts 2 (q !6 -|- 0!y L Oio) 
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Again we can use the relations derived from the C, P and T to relate the helicity 
amplitudes and thus to check our results. These relations are: 

\Mx 1A2A3A4 I |-^A 2 AiA 4 A 3 | 

p 

|dT^jA2A3A4| Ai —A2 —A3 —A4 I 

T 

|.^AiA2A3A4| |-^A3A4AiA2l 

In order to relate some helicity amplitudes we have also taken into account whether 
the scattering particles are identical or not., There are still 15 independent amplitudes 
remaining after applying these symmetry relations. 

The W~^W~ reaction represents another source of Z^Z^ pairs from pp beams. 

Every helicity amplitude in this process receives contributions from the t and u channels, 
by exchanging a.W or a vr, and directly from the vertex W^W~ Z^ Z^ (Figure Therefore 
we can write: 


^AiA 2A3A4 — + Mt2 + + ^U2 + ^ 4 ] AiA 2 A 3 A 4 

where the subscript 1 ( 2 ) is referred to the exchange of a gauge boson (Goldstone boson) 
through the corresponding channel. There are 81 helicity amplitudes, but the number of 
independent MA 1 A 2 A 3 A 4 is reduced to 25 by means of symmetry relations derived from P 
and C invariance. 

As we expected, the amplitudes that we have obtained satisfy these equalities. On the 
other hand, all the chiral parameters, (but aio), affect this process as we can deduce 
from the Feynman rules eqs.(^ and ^) and from the different diagrams that contribute 
to the collision W+W- Z°Z°. 

3.2 final state 

As it was previously said, we have also studied the W^Z^ hnal state. The different sources 
for W^Z^ pairs in pp colliders at tree level order are the following: 


qq' 


qq 


W^Z^ 


w^zl 


w^z^ 


wtzl 
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W^Zl W-Z^ 

W£Z° W-Z° 

WEZl W-Z° 
w^7 ^ W+Z° 
iy+7 ^ w+z° 

W-J W-Z° 

WE7 w-z° 

where we call q = u,c and q' = d, s, b. 

Most of the boson pairs W^Z^ are obtained in pp colliders via quark-antiquark an¬ 
nihilations. Using the chiral Lagrangian framework we calculate the corresponding cross 
section at tree level order. As we can see in Figure H, three standard s, t and u diagrams 
contribute, but the new physics coming from the non-linear lagrangian is isolated in the 
three boson vertex, which in this process affects only the s channel. 


M = Ms + Mt + Mu 


If we look at the Feynman rules (Figure]^) we can deduce that the set of ak parameters 
contributing to this process are ai, 02 , 03 , ag, Og, an, ai 2 and ais eq. (PT|) . In the 
hadronic collider case, the three boson vertex has been studied in this process with chiral 
Lagrangians, in at tree level, using the ET (by J.Bagger, et.al. in |^), or in our 
previous study [|T^ for the case g' = 0 and including the running of the couplings. In 
this last case we analyzed the sensitivity to the ag parameter. Here we work hrstly at 
tree level order, but including the dependence on the ak parameters in all the channels 
to produce W^Z^ pairs in pp colliders. Finally we will include the running parameters 
effect in section 8 . 

Other mechanism to obtain W^Z^ pairs is through W^Z^ or collisions transver- 
sally or longitudinally polarized at the initial state. In both cases the different polarization 
amplitudes are obtained by adding the contribution of the s and u channels and the four 
gauge boson vertex in (Figures | and |T^. The exchanged particles in the s and u diagrams 
are a gauge boson (si, uE) and a GB (s 2 , ^ 2 ), respectively. 


AfAlA2A3A4 — [^Sl + ^32 + + Mu2 + M 4 ] A 1 A 2 A 3 A 4 

In the W^Z^ W^Z^ collisions there are initially 2x3^ helicity amplitudes. This 
2 factor can be dropped using C invariance that relates the magnitude of the helicity 
amplitudes Mx^x^x^x^ in the process with hF+ and with W~. If we apply the relations 
between amplitudes derived from P and T invariance, there will be only 25 remaining 
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amplitudes. These subprocesses cross sections are affected by all the but aio, as we 
can deduce from the Feynman rules eqs.(|^ and |3^ . 

The number of helicity amplitudes corresponding to the processes hF^y — 
is 2 X (2 X 3^). The 2 factor disappears when we take into account C invariance as in 
the previous case. The number of remaining amplitudes is reduced again, by another 
2 factor after applying P invariance. Therefore in these processes we have 3^ = 27 
independent helicity amplitudes that we calculate using the Feynman rules shown in 
eqs.(^ and Here we can observe that the set of the chiral lagrangian parameters 
affecting the fF^y —cross sections is the following: oq, Oi, 02 , 03 , as, ag, an, 


ai 2 cind a^s. 

In order to see the relative importance of the different channels we have evaluated the 
contribution to the total cross section to obtain Z^Z^ and W^Z^ pairs in pp colliders 
from the different considered subprocesses in a typical example. We have chosen the a^ 
parameters that mimic the MSM with a heavy scalar Higgs whose mass is tuh = ITeF . 
The rates for the contributions of the different subprocesses to the total Z^Z^ final state, 
qq, gg, and W^W~ are, respectively: 45, 17, 3, 35% and to the total W^Z^ the 

rates are 61, 22, 17% coming from qq', W^Z^ and W^y fusion. 

As it has been discussed above, the essential point in our calculations is that we have 
not used the ET P, M to obtain the scattering amplitudes of the gauge bosons as it 


is customary. We have calculated the corresponding tree level amplitudes for the gauge 
bosons and projected them into their transversal or longitudinal components. Therefore 
we do not have the limitations that appear when the ET and the chiral lagrangian for¬ 
mulation of the ESB are used together IQ. The chiral lagrangian approach provides a 
low-energy description of the GB dynamics as an expansion on the momenta over Attv and 
the ET refers to the large energy relation between GB and the longitudinally polarized 
gauge bosons S matrix elements. In our case we only need to fix an upper energy bound 
(as it was said i^max = 1.5TeF) so that we can safely apply yPT. 

However, the ET can be used as a helpful tool to test our calculated longitudinally 
polarized amplitudes. In order to check our results we have compared our tree level 
expressions with those calculated applying the ET to the corresponding Goldstone boson 
amplitudes. We have found both amplitudes to agree at lowest order in g and g' for the 
cases Z^Z^ —>• Z^Z^, kF/lFh —>• and W^Z^ —>• hF^Z° as it was expected according 

to the results of [10|. However, the details of this comparison will be described in detail 


elsewhere since they concern the applicability of the ET which has not been used in 
the present computation. Here we only wanted to quote that our results are compatible 
with the information that one could have obtained from the ET on the above mentioned 


processes. 
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4 Signatures and background 

As it was stressed in the previous section we have calculated the tree level scattering 
amplitudes, up to order corresponding to all the helicity states of the gauge bosons, 
without using the ET. This fact allows us to include the dependence on the ak parameters 
in all the helicity amplitudes. In contrast, when the ET is used, the only channel where 
the ak coefficients are taken into account is the following: 

VlVl^VlVl (V = Z°,W^) 

The dehnition of this channel as the signal was considered appropriate because the 
VlVl —> VlVl channels are expected to be strongly interacting at high energies, if the 
GB are, due to the ET. However, our computation makes it possible for the hrst time to 
study the effect of a strongly interacting ESB in other gauge boson polarization states by 
means of the ak parameters. In fact, the aim of this work is to see how measurable will 
these parameters be at the LHC. In order to make this point more precise we must dehne 
the statistical significance corresponding to some given value of the chiral parameters ak- 
With this goal in mind we dehne our signal and our background as: 


(Signal) ns = |iV({a}) - iV({a°})| (10) 

(Background) ub = iV({Q;°}) (11) 

where iV{Q;}, are the total number of or W^Z^ pairs obtained for some 

given experimental cuts when the chiral parameters have been set to the values {a} or 
{q;°}. The background is dehned in terms of some reference model For simplicity 

this model has been taken as the one with all the parameters set to zero i.e. a® = 0. We 
call it Zero Model and, incidentally, it corresponds to the MSM with an inhnite Higgs 
mass. 

For the hnal state Z^Z^ we have considered the following processes that contribute to 
the signal and the background: 




In addition we include the Z^Z^ production via gluon fusion and qq annihilation which 
do not depend on ak and therefore contribute only to the background. The experimental 
signature for this process consists of four leptons as the result of the Z^Z^ decays. This 
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gives a clean and distinct signal because the pairs can be fully reconstructed. The 

disadvantage is the rather small leptonic branching ratio (0.44%). 

In the W^Z^ case we have considered the channels: 

W^Z^ W^Z° 

gg' 

We focus our attention on the gold-plated events where the IV^ and Z^ decay to the 
charged leptonic hnal states (/ = e,/i, z/). The corresponding branching ratio is 1.3%. 
All these cross sections depend on the chiral parameters so that they are taken into 
account in the signature and background calculations. 

The main source of Z^Z^ or W^Z^ pairs in pp colliders is via quark-antiquark annihi¬ 
lation. The total production rate of gg —>■ Z^Z^, in the studied cases, is 20—50% than that 
from gg Z^Z^ depending on the top quark mass (we have chosen m* = 170 GeV). On 
the other hand, the production rates of Z^Z^ or W^Z^, via gauge boson fusion are sup¬ 
pressed by powers of (a/sin^6'^) due to the application of the W effective approximation 
|T5| , pTj to obtain the initial bosons from pp beams, as we will see in next section. 


5 The proton and gauge boson structure functions 

Here we describe how to compute the total cross sections of the different processes studied 
to obtain Z^Z^ or W^Z^ pairs in pp colliders. 


pp {gg V 3 V 4 ) X 

pp ^ {gg^Z^Z^)+X 

PP ^ {V,V2^VsV,)+X 

We have to integrate the differential cross section for the subprocess, {da/dcos9), with 
the distribution functions of the quark, antiquark and gluon (given by /*, fj and g) inside 
the proton: 


a{pp {gg' V 3 V 4 ) + X) 


J J J dxidx2dcos9fi{xi,Q‘^)fj{x2,Q‘^) 


dd 

dcos9 


{qq' ^ % 3 % 4 ) 
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a{pp ^ {gg ^ Z^Z^) + X) — J J J dxidx2dcos9g{xi,Q‘^)g{x2,Q^) 

-^(gg Z«Z») (12) 

acosU 

These formulae are used to compute the processes: 


W^Z^ 

In order to compute the number of events of Z^Z^ and W^Z^ produced in pp collisions 
via gauge boson fusion we apply the effective W approximation and we use the 

formula: 



a{pp^{y^^2^V^V^)+X) 


J J dxidx2dcos9fi{xi,Q‘^)fj{x2,Q‘^) 



d^L dd 
dfdfj dcos9 






(13) 


Thus the total cross section in pp colliders (like LHC) can be written as the result of 
the convolution of the subprocess cross section with the V 1 V 2 pair luminosity in pp beams. 
This luminosity is calculated from the convolution of the double bremsstrahlung of the 
(V 1 V 2 ) from the quark structure functions. Thus, d^L/dfdfj is the luminosity function 
for the gauge boson pair V 1 V 2 to be radiated from the quark pair qiQj. It depends on 
the helicity state, transversal or longitudinal, of the initial bosons Vi and V 2 . Therefore 
we have to separate in our computations the contribution of the different polarization 
channels of the initial bosons. 




W+W- 


w^z° 

^ w^z^ 

-H 

w^z° 


The amplitudes and differential cross sections for all these processes have been obtained 
as it was described in the previous sections. According to the effective W approach we 
take the following functions corresponding to the probability of a gauge boson Vr or Vl 
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(V = W^,Z^) or a photon, to be radiated from the quark q, with a momentum fraction 
of the quark, x. 


fq/VT{x 

fg/V^ix 

fq/'ri.X 


^ x^ + 2 {l-x)^ ( 
fv -^-In 


2 x 


K, 


) = fv- 


X 


X 


^, l + (l-x)^ ln IE 

2tt ^ X \m 


V = w,z 


(14) 


where the values of the fy.'g depend on the particular gauge boson as well as on the 
type of quark that it comes from 


fw{x) = 
fz^^{.x) = 

fzAx) = 


a 


At^x 


w 


a 


167ra;iy(l — xw) 
a 

167ra;pv'(l — xw) 


1 + (1 - 

1 + (1 — -XwY 


where xw = sin^6*w. 

The variables r and rj are related to the momentum fractions of the quarks by 
3 ^1,2 = The connection between the variables f, fj and the momentum frac¬ 
tions of Vi, V 2 respect to gi, Qj, xi and £2 is given by Xy 2 = Note that the rate 

of transversally polarized gauge bosons obtained is enhanced by the logarithmic factors 
eq. (p^ with respect to the longitudinal gauge boson production. 

The structure functions of quarks and gluons we use are those of EHLQ set II 
with A = 290MeV and we neglect the contribution of the top quark to the proton sea. 
However we have also studied the effect of changing the structure function on our results 
(see below). The assignment for the appearing in the distribution functions is = s 
for qq' and gg processes, and for elastic gauge bosons scattering processes. The 

resulting integrals are computed using the VEGAS Monte Carlo program . 

With the described machinery we have built up a big code which computes the total 
number of expected gold-plated events (for some given LHC integrated luminosity) in 
terms of the chiral parameters and the hnal state cuts . 

Initially we impose the following cuts on the invariant mass of the weak boson pair, 
the transverse momentum of the hnal Z^ boson and the rapidities of both bosons: 
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200GeV < Vs < l500GeV 

Ptz > IQGeV 

\yimJ = \y2^J = 2.5 

From now on, we will call them the minimal cuts. The upper limit in the invariant mass 
ensures the validity of yPT approach. For the rest of the paper we compute the number of 
and W^Z^ obtained in the LHG with an integrated luminosity of 3 x This 

approximately corresponds to a total LHG working time of 1 year (3 x lO^sec) assuming 
a luminosity C = 


6 Parameter sensitivity, optimal cuts and measur¬ 
able parameters 

In previous sections we have presented in detail how we have computed the number of 
Z^Z^ or W^Z^ pairs obtained at the LHG. If we fix the integrated luminosity, the pp 
center of mass frame energy for the LHG {y/s = 16TeV), and the upper bound on the 
invariant mass at the subprocesses (V^max = l.STeld), the total cross sections will depend 
on the chiral parameters and on some kinematical cuts. For the sake of simplicity we 
will only study the subset of the coefficients which is needed to reabsorb the one-loop 
divergences obtained calculating with the chiral Lagrangian to lowest order. They are Oq, 
ai, 0 ^ 2 , as, a 4 and as ||^ (however, our code includes as well the contribution of the other 
chiral parameters). There is also the possibility to modify our results choosing different 
kinematical cuts, q (ci = Vsmin, C 2 = Przmin 


and cs = = y2^V- Thus we can define 


two vectors: a and c, as follows: 


a — (ao, ai, a 2 , as, a4, as) 

C = (Ci,C2,Cs) (15) 

which are inputs of our code; whereas the number of Z^Z^ [i = 1) or W^Z^ [i = 2) 
events is the output. Therefore we write the result of our computations as N^^\a;c). In 
principle one could expand these functions around the Zero Model (we defined it previously 
as {a° = 0,V/c}) and the minimal cuts defined above so that one could write: 


V(*)(a;c) = V(')(a°;c““)+ ^ 


k=0 


dN^i) 

dai 


'“fe 


=o«fc + 


1=1 


dci 


^.(Q-cr) 


\Ci=C 


+0{V) + 0{{c-c'^V ) 


(16) 
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If the dependence of the nnmber of events on the parameters and the cnts were ap¬ 
proximately linear, this formula could be used to compute N^^\a;c). However, we will 
see that this is not always the case. Indeed, c) is a polynomial in and nonlinear 

terms can become important even for moderate values of these parameters. 

In order to see how the number of events changes when one of the parameters varies 
by some amount Aa^, we £x the kinematical bounds q to the minimal cuts and set the 
other chiral parameters to zero. We choose the following set of values for Aa^: 

Attfc = ±10“^ ±5 X 10“^ ±10“^ 

The result of our computations can be found in Figures 0, to 0. We can also dehne 
the sensitivity function (c) associated to the ak parameter, with the kinematical cuts 
c, and for the hnal state i, as follows: 


’?(c) = 


awh) I 


a; c) 


dai 


lofc=0 


(17) 


The different sensitivity functions for the minimal cuts that we have obtained are 
displayed in Table Obviously, the |s^*^(c)| values are a measure of the variations of 
the number of events with ak- However, they are not a direct measure of the statistical 
signihcance of the corresponding parameter variation. 

We can also study the number of and W^Z^ event distributions with respect to 
the cuts given in the minimal invariant mass (ci), the minimal transverse momentum of 
Z^ in the hnal state (02) and in the maximal rapidity of the hnal bosons (03). In these 
computations we have hxed the a parameters at their values in the Zero Model, and the 
other two kinematical cuts have been set to their minimal values. The results obtained 
in this way are shown in Figures 0 to |T^. 

In order to perform a statistical analysis, we dehne the function as: 




c = 


_ |wW(afe;c)-ArW(a°;c)| ^ ns 
^N0){o°\c) 


(18) 


The function is a measure of the statistical signihcance of the signal (corresponding 
to the increments Aa^) relative to the background, provided N^^^ak', c) is large enough to 
apply the Central Limit Theorem. In this case is an estimate of the number of sigmas, 
and therefore it dehnes the conhdence level for the hypothesis that ak is diherent from 
zero. In order to evaluate we use the number of events that was previously obtained. 
The resulting functions are shown in Figures to 

From Figures 0 and 1 ^ it is clear that there is no linear dependence of N^^^a; c) on 
Cl and C 2 . Moreover, looking at Figures 0 at 1^ and from Figures ^ at |^, nonlinear 
behaviour is observed with respect to ai, 02 and 03 (for W^Z^ hnal state) nor when | 04 | 
0)5 I are bigger than 0.001 in both channels, Z^Z^ and W^Z^. 


or 
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By looking at Figures]^ to we can see that the slope corresponding to the different 
cnrves varies according to the different valnes of and the channel 

or W^Z^). Indeed, the greater is the slope, the higher is the statistical significance. 
Moreover it is already possible to observe that some of the a parameters have no chances 
to be measnred at the LHC, at least in the way described here. For this reason, in the 
following we will concentrate our statistical analysis on the set of a parameters which 
can be considered as potentially measurable. Being more precise, we define a parameter 
ttfc as potentially measurable in a given channel when the corresponding to the 

maximal considered variation (| |= 0.01), is bigger than 0.5. According to this we 

find that the potentially measurable ak are the following: 

• 0:3 < 0 and 0:3 > 0 (both in Z ^ Z ^ and W ^ Z ^ channels) 

• 04 < 0 (IF^ZO) and 04 > 0 (Z^Z^) 

• as < 0 (PF^ZO) and as > 0 {Z°Z^) 

In the following we carry out an statistical analysis with all these potentially measurable 
parameters in the corresponding final state. To do so, the next step is to look for the 
optimal cuts in the minimal invariant mass and transverse momentum for the detection 
of a certain Aa^ = (a^ — a°) different from zero. Initially we impose the minimal cuts on 
the kinematical variables. The optimization procedure we use is to search for the pair of 
cuts (\/^, Ptz) kinematical region 

200GeI/ < \/^ < IbOOGeV 
lOGeV < ptz < \J Smax/4 - ml (19) 

that maximize the function r^\ak', \/^,Ptz) defined as follows: 


r^k\ak;^\pTz) = 


|AW(afc; Vs > y/¥,pTz > Ptz) “ Vs > y/¥,PTZ > Ptz)\ 


Ah) I 


a 


s > y s^,pt > Pt) 


( 20 ) 

To find the optimal cuts for the different parameters and channels we build a bidi- 
mensional grid on the plane (\/§,Ptz) so that the different points are separated from each 
other by the increments (Aa/^ = 50GeF, Apx = 50GeF). For each pair of points in this 
grid ( a /^, Ptz) compute the total number of events in the Zero Model and in a ESB 
scenario corresponding to a certain positive or negative Aa^. Finally, we evaluate the 
rl^\ak', \/^,Ptz) function. The optimal cuts, c°p, are those which maximize eq. (|20|) , and 
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therefore, the conhdence level. Of course, these c°p which have been found in this way de¬ 
pend slightly on the particular choice of Aa^ for each parameter. In our computation 
we have taken the typical values: 


Attk = ak — = ak = ±0.005 (21) 

The results are collected in Tables ^ and There it is displayed, for each Aa^ in the 
corresponding channel, the optimal cuts, the number of events that satishes these cuts 
and the statistical signihcance function obtained with the minimal and the optimal 
cuts. In all cases we have hxed a maximal rapidity of 2.5. 

The c°P we have obtained by maximizing the function eq. (l20|) will be considered 
from now on as the optimal cuts to detect either a Aak < 0 or a Aak > 0. For example, 
if we wanted to observe the signature corresponding to a heavy Higgs S'M-like ESB 
scenario, the best choice would be the hnal state with the optimal cuts obtained 

for Acts > 0 (A/smin = 1150GeH, PTZmin = 400GeH). This is so because, at the tree level 
order, the only ak different from zero needed to mimic this scenario is cts (as = (8m|^)). 

As it can be seen in Table the optimization procedure has clearly improved the 
statistical signihcance in most of the studied cases. Only when we tried to optimize 
the signature corresponding to a as > 0 and as < 0, both in the hnal state, we 

obtained no signihcant improvement with respect to the minimal cuts. Therefore we hnd 
it impossible to detect a signature of a negative or positive as in Z^Z^ events. 

Finally, we have carried out a similar optimization procedure for and as in which 
we searched for a slightly diherent pair of optimal cuts: (^Smax, PTZmin)- We obtained 
c°P = (1500GeI/, 500GeV) for = -0.005 (fF±Z°), = 0.005 {Z^Z^) and as = 0.005 

[Z^Z^) and the function was respectively 6.02, 1.56 and 3.82. On the other hand, 
we also found the optimal cuts to detect as = —0.005 in the hnal state (c°p = 

(1500GeI/, 400GeH)) and the reached statistical signihcance was 2.97. These results are 
at most equal to those of Table 

Moreover, we have repeated our whole study to hnd other sets of optimal cuts belong¬ 
ing to the parameter space, such as (V^max, PTZmin)^ (V^max, PrZmax), (V%^, PTZmax)- 
In all these cases, the statistical signihcances obtained are smaller than those correspond¬ 
ing to the hrst choice of c°p (c°p = (V^mm,PTZmin))- Therefore, in the following we will 
apply this type of optimal cuts. 

As it was mentioned in the introduction the main aim of this work is to determine 
which chiral Lagrangian parameters will be more easily measurable at the LHC and the 
size of their corresponding statistical errors. Initially, we made a criterion to choose a 
set of potentially measurable ak parameters, on which we have concentrated our analysis. 
After carrying out the described optimization procedure, we conclude that the statistical 
signihcance to measure most of these chiral Lagrangian coefficients, has clearly increased. 
In the following we will treat only these measurable ak- The set of the ak parameters 
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obtained as the measurable ones are contained in Tables ^ and Here we exclude the 
study of tts in the hnal state, as it was previously argued. 

In order to obtain an estimation of the statistical errors we will work under 

some hypothesis. We will assume a linear behaviour in ak and that the number of events 
is large enough to consider following a gaussian distribution. 

At this point, we dehne the statistical error associated to a certain as 

the minimal value of this parameter that could be detected in the LHC with a statistical 
signihcance, eq. fpAf ), equal to one sigma, after applying its corresponding optimal cuts. 
Using the results displayed in Table ^ which satisfy the optimal cuts, we obtain the 
(Aafc)g^g^^ in this way, and we show the results in Table As it could be expected, 
the smallest (Aafc)g^ 3 ^^ value corresponds to the biggest statistical signihcance, which is 
(Aa 4 )g^.g^^ = 8.31 X 10“"^, in the hnal state. 


7 The effect of the structure function indetermina¬ 
tion 


In all the previous calculations we have used the EHLQ structure functions (set 
II with A = 290MeU) to obtain the total number of and VC^Z^ events. Here we 
wonder how our results could be ahected by the indetermination in the parton distribution 
functions. In order to estimate this ehect we have selected two new parametrizations of the 
structure function: the set MRSD— (Martin, Roberts and Stirling |^) and the GRVHO 
(Gliick, Reya and Vogt [^). Both them exhibit a similar behaviour of sea quarks and 
gluon distributions at low x, that grow when a; ^ 0. The values of A considered in 
these two sets are the following: Amrs = 230MeU, Agrv = 200MeU. With these new 
structure functions we evaluate the total number of Z^Z^ or W^Z^ events, in the Zero 
Model with the different optimal cuts. The way we proceed is to compare our estimation 
of the number of standard deviations obtained by changing these parton density function 
sets (we will call it with the statistical signihcance corresponding to a certain 

signal, after having applied the optimal cuts (the r^k\c°^) functions shown in Table 


We dehne rgJruc. similarly as in eq. (pAD : 


Ah 

struc 





( 22 ) 


Now, to obtain we need to evaluate the total number of events in the Zero Model, 
with certain kinematical cuts, c, but using two diherent parametrizations called set and 
seUfor the structure functions. 


23 








In Table ^ we show the results of our calculations. We represent the total number of 
and W^Z^ events obtained for the optimal cuts, when the EHLQ (set II), MRSD — 
and GRVHO sets are used, as well as the number of sigmas Ustruc. eq.(p2D corresponding 
to change one of these sets with respect to the others: 

• Tstruc. = ’"12 {EHLQ with respect to MRSD—) 

• Tstruc. = ’"13 {EHLQ with respect to GRVHO) 

• Tstruc. = ’"23 {MRSD— with respect to GRVHO). 

Now, we pay attention to the r*, functions in Table ^ and Ustruc. in Table ^ and we 
compair them. We will say that the typical chosen value of ak {c^k = —0.005 or = 0.005) 
can be measured in the LHO with a statistical signihcance given by r^, if the following 
inequality is verihed: 


’"fc ^ ^ ’"struc. (^ 0 ) 

When we make Ustruc. = ’"12 or Ustruc. = ’"13, we must reject the value as = 0.005 in 
Z^Z^ and W^Z^ channels, and as = —0.005 in hnal state, because relation (PB| ) is 
violated. However, if we take rstmc. = ’"23 we can measure positive and negative values of 
as, ctA and as parameters in some Z^Z^ or W^Z^ hnal state, with a statistical signihcance 
given by r^. On the opposite, for as = ±0.005 in Z^Z^ hnal state, eq.(p^ is not fulhlled, 
so we can only measure these values of as through the W^Z^ channel. This last choice 
(’"struc. = ’"23) is most appropriate since MRSD— and GRVHO agree much better with 
some recent experimental results p6|] obtained at HERA. In any case, it can be expected 
that future new experimental data coming from HERA and even from the LHO itself, 
will allow US to reduce the size of these uncertainties Tstruc.- 


8 Renormalization effects 

The other ehect we are going to take into account is the dependence of the ak parameters 
and the couplings g and g', on the energy. Our whole previous study was carried out 
at tree level order. Now we want to include some quantum ehects by means of the 
Renormalization Group Equations {RGE). By that we will understand the one loop 
contributions coming from and /Iym plus the tree level contributions coming from 
C[ ± to the corresponding beta functions. 

As it is well known from the RGE, the renormalized Green function with the renor¬ 
malized parameters A and m at the renormalization scale /i and the same Green function 
with the parameters at another scale pe*, are related by the equation: 
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Gk')?!,A. ji) = GS5‘’(pj,A(t).(ie')exp[-n f i{\(t'))dtf] 

Jo 

Where we include in A all the ak parameters and the electroweak coupling constants, 
g and g'. Thus the renormalized parameters are scale dependent. The A stands for all the 
running coupling constants that depend on the renormalization scale. In order to obtain 
these A we need to calculate the and 7 (A) functions, following the well known RGE 
techniques: 


m 

7(A(t)) 


dt 

1 dlogZsv 

2 dt 


(25) 


At present, when using the effective theory to describe the ESB sector we only know 
these functions to one loop-order. The exponential factor in eq.(P^) is the anomalous 
dimension term, that depends on the wave function renormalization. We have estimated 
the order of magnitude of this effect to calculate the statistical signihcance r^, in a typical 
case, with a subprocess energy of ITeV. We obtained changes on of 0.1%. Therefore 
they are completely irrelevant. Thus we neglect this contribution and we will only consider 
the dependence on the renormalization scale, p, of ak, g and g', in eq.(^^. We refer to 
this approach as the tree level approximation improved by the RGE. 

Now, we have to calculate the running coupling constants ak{g)-i and g\g), and 
replace their tree level values in our cross section formulae by the corresponding energy 
dependent running coupling constants evaluated at the subprocess center of mass energy 
(we have taken /i = a/^). 

In order to do so we begin with the ak parameters and we use the dimensional regular¬ 
ization scheme since it is the most appropriate for gauge theories as well as for non-linear 
sigma models. In []^, Longhitano obtained all the divergences that appear in one-loop di¬ 
agrams calculated with the effective lagrangian at lowest order + Cym)- They could 
be absorbed by redehnitions of a subset of the ak parameters {k = 0,1, 2, 3,4, 5). By 
means of the dependence of ak on the e parameter [|^ 
ak running expressions: 


it is easy to obtain the following 


a^(p) = a^(po) - Y^^tgXlog 

af(/i) = af(po) - Y^^tg^wlog 
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= a^ino) - 


= a^(/io) 


Rt 


levr^ 12 
1 1 




16 v? 12 ^°S (7. 


af 7) = af 7o) + 


(5) 


af7) = «f7o) + TTZ2l^log 


/ ^ 

167r2 12*"^ /i ^ 


( 26 ) 


Apart from these equations, we need also to calculate the running of the gauge cou¬ 
plings ( 7 i?(/i) and 7(/^)- usual we have 


9R = ^^0 

/ /7'/1/2 / 

S'r ~ ^ 3 % 

The only differences with respect to the linear model, given by 6Z3 and 6Z'^ functions, 
come from the scalar sector contribution. In our case, they are 


<5^3 


9 


IGtt^ 6e 

V2 1 


9 


167r^ 6e 


Taking into account loops of scalars, fermions, ghosts and gauge bosons coming from 
the rest of the SM particles we arrive to the following values of the (3 functions eq. (P^D , 
f3g and (3gi: 


where 


Pg = -Cg9^ 

7 = -Cg'9" (27) 


Cg 

Cg' 


1 13 

IGvr^ 4 
1 27 

IGvr^ 4 
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The eqs.(^) are standard corresponding for gange theories. The hrst one leads to 
asymptotic freedom since Cg > 0, whereas the second corresponds to an abelian gange 
theory. Therefore one has: 


= 


= 


9Uf^o 


1 + 


1 13 

167r2 2 

,/2 


4(/io)log^ 


9 ' 


1 - 


1 27 „,2 

Ww^ 2 


9'K(Ao)log 


(28) 


Mo 


Once we have obtained the expression giving the dependence of the ak, g and g' on 


the renormalization scale eq.(^), we can stndy cnantitatively this effect. The way we 
proceed is the following: 

- First we factorize the g and g' conpling constants so that the cross sections appear 
with the same power in the electroweak conplings than in the MSM. We are 
referring only to the dominant terms. Moreover, as we can see in the Appendix, 
eqs.(^, ^ and ^2]), higher powers on these conpling constants appear as factors of 
the ttfc parameters. 

- In all onr calcnlations, for the sake of simplicity, we are considering the physical 


masses 


27|: mw = 80.6GeId, mz = Ql.lGeV and = 0 (except to in gg —>• 


process where we have taken a top qnark in the loop with a mass of 170GeI/). 

We replace in the cross section the constant valnes oi ak, g and g' by those which 
are renormalized and scale dependent, a^(/ii), gniju) and g'fi(/u) given by eqs.(|26| and 


). We take /i = ^/S, where is the invariant mass (the center of mass energy in 
the considered snbprocess) and po = mz- Thus we have substituted the tree level 
values oi ak, g and g' by the running coupling constants depending on the energy 
scale of the subprocess. 

The assignment for the numerical values of gji(/uo), g'ni^Ho) and akif^o) parameters 
appearing in our formulas has been done as follows: 

- The quantities g^ and g'j^, at the chosen scale /xo = have been taken from recent 


LEP data 27 


- The symmetry breaking pattern that we have considered corresponds to neglecting 
all the af parameters at the /xq scale. We refer to this scenario as the Running Zero 
Model'. 


akiRo) = 0 (A: = 0,1,2, 3,4, 5) 

ak =0 (A: = 6,7,8,9,10,11,12,13) 
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After all these changes in the program, we apply the different sets of optimal cuts 
previously obtained, (Table 3) and calculate the total number of events to produce 
or W^Z^. The corresponding results are displayed in Table |^. We present for each pair 
of optimal cuts, the total number of Z^Z^ or W^Z^ events, in the Zero Model and with 
g and g' constant {N^''l{a^]g,g',c)), and in the Running Zero Model with gniVS) and 
9k(VI) (/V«( a^{\/§); gii{'/§), g'j^{\/§), c)). We also give an estimation of the statistical 
signihcance or the number of sigmas, , corresponding to include the dependence on 
the energy of the parameters, with respect to the tree level results. Here, we calculate 
rj^^^ function analogously to eqs.(|^ and (^) 



lN(^)(a°;g,g',c) - N(^^(a^(\/§);g(\/3),g'(Vs),c)l 
^/N(^)(aO;g,g',c) 


(29) 


As it can be observed in Table |^, the function, that reflects the relevance of this 
effect, is very high when we use the minimal cuts. Besides, with all the optimal cuts, 
the statistical signihcance function whose values go from 1 to 4 standard deviations 
is also important. Thus, the effect of the running of the couplings is important and it 
should be taken into account if one wants to compute the total number of Z'^Z^ or W^Z^ 
events. 

Now one could ask how the r^ that were previously obtained (Table ^ would change 
if we used the improved tree level approximation with the running parameters, instead 
of working only at tree level order. In sight of the results contained in Table and if 
we supposed linear behaviour we would obtain huctuations in the statistical signihcance 
functions, r^, varying between 2 to 15% with respect to their values calculated without 
including the dependence on the energy of the parameters. 

We can summarize these results saying that the inclusion of the running coupling 
constants in our calculations leads to important variations in the total number of events, 
given by functions. However, this ehect is not so relevant when we want to determine 
the Tfc functions corresponding to each since the diherences obtained with respect to 
the tree level results are small. Therefore, the running of the couplings does not change 
signihcantly our previous discussion about which couplings will be measurable at the 
LHC. 


9 Conclusions 

Using yPT to describe the ESB sector, we have studied the different processes contribut¬ 
ing to Z^Z^ and W^Z^ hnal states at the LHC, considering only gold-plated events. The 
main novelty in our analysis is that we do not use the ET but we have worked explicitly 
with all the the gauge boson polarization states. Thus we can study the low energy region. 
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after having imposed a maximal bound on the subprocess energy of l.hTeV. We have 
elaborated a FORTRAN code that generates the subprocess cross sections at tree level 
order, producing or W^Z^ pairs. 

With this code we have carried out a systematical study of the possibilities for mea¬ 
suring the chiral parameters at the LHC. This analysis includes an optimization 
procedure, that makes possible to obtain the greatest statistical signihcance for mea¬ 
suring certain values of these parameters. We have also computed their corresponding 
statistical errors, (Aafc)stat .5 set the minimum values of these to be unambiguously 
detected with a certain statistical signihcance at the LHC. The results of our studies are 
displayed in Tables]^ and From our results it is clear that the Oq, Oi and 02 parameters 
cannot be probed at the LHC. On the contrary, as it can be seen in Tables ^ and |^, 
the Tfc function obtained for 04 and can be, for some channels, greater than 3 sigmas 
(2 sigmas for a^). In fact, the parameters that could be more easily measured in the 
LHC are 04 < 0 in W^Z^ and as > 0 in Z^Z^ channel. As it can be expected their 
corresponding statistical errors are the smallest. In all our calculations, we have hxed a 
running time corresponding to one full LHC year according to the nominal luminosity 
L = (Tint. = 3x 10^pb~^). The results corresponding to different integrated 

luminosities can be obtained just rescaling the statistical signihcance by the square root 
of the running time. 

On the other hand, we have also estimated the size of the imprecisions in our calcu¬ 
lations due to the indetermination in the structure functions by choosing three diherent 
distribution functions. The results (which can be found in Table ^ show that the errors 
coming from this ehect are smaller than the statistical errors obtained for the a parame¬ 
ters. Thus our present ignorance about the parton distribution functions does not restrict 
the measurement of the chiral parameters. 

We have also tried to improve our tree level results considering the dependence on 
the subprocess energy of the weak couplings and the chiral parameters. The size of these 
effects is shown in Table From our results it is clear that the total number of events 
changes when this effect is taken into account. However, our previous conclusions about 
which parameters can be probed at the LHC, do not change at all. 

Moreover one could ask about the possibility of increasing the number of measurable 
chiral parameters by relaxing some of the experimental assumptions of this work. For 
instance we could also take into account the gauge boson hadronic decays, separate the 
hnal polarizations, etc... In that case we could certainly enhance the statistical signihcance 
and more chiral parameters could be probed. However it is not possible today to have 
any idea about how well this new information will be obtained at the LHC and for this 
reason we have not considered such possibilities. 

Finally we have also assumed that the or W^Z^ pairs produced via electroweak 
gauge boson fusion cannot be separated experimentally from those coming from another 
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sources such as gg fusion (in case) or quark-antiquark annihilation. Nevertheless 

some forward calorimeters could, presumably, be incorporated to the LHC detectors. This 
fact could allow for a jet tagging at a certain level of efficiency. This experimental im¬ 
provement is a realistic way to increase the statistical signihcance functions to measure 
the chiral Lagrangian parameters at LHC. Work is in progress to treat this possibility 
and to analyze how much the sensibility to the parameters could be enhanced. 
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Appendix 


In this Appendix, we write some of the Feynman rules obtained from the quantized 
chiral electroweak Lagrangian £nlsm eq.(0), described in section 2. 

The different contributions to £nlsm are given in eqs. flai and IP , where the Landau 
gauge has been chosen. 

We show the Feynman rules corresponding to the propagators and vertices used in our 
calculations: 

i) The TT and propagators in the Landau gauge have the following expressions 
(see Figure |^): 


-zA(7r±) = ^ 

ii) In the following, we write the Feynman rules corresponding to the vertices 
with either 3 gauge electroweak bosons (IviVaVs) or to 2 gauge bosons and a GB (IviVjtt) 
displayed in Figure 


G+w-zo 




^w-w+z° 


^W^W+Ax 


—sec9wg^,u{sm'^0w - + — [{q■r)g^,u - q^,r^]{aism9w - a8Cos9w + 

^ X) 

2g‘^ 2g‘^ 

ai3sec9w) + —a3sec9w[{k-r)g^,y - k^r^] + — [{k.q)gf,„ - qf,k^]{a2sm9w - 
2g‘^ 2g‘^ 

(tta + ag)cos9w) H- aiiCos9w[kij,qu - r^j^kiy] - allsm‘^9wsec9w[k‘^g^_lu - 

g^ 2 ^^ 

k^,qu + r^r^] + ^al2sec9w[k‘^g^,u - k^q^ + + q^q^] - -^ai3sec9wqf,qu 

^2 


g V . 2 
-sm 


2 g‘^ 2g‘^ 

Owgfiu - [{q■r)g^,v - q^r^]{aiCOs9w + agsin^w)- [{k.q)g^y 

^ V V 

2 g‘^ 

-qf,ku]{a2Cos9w + (as + a9)sin6'ty) + -^aiism9w[k‘^giMu - k^r^ + 
igcos9w{l + a3g\ec‘^9w)[{k - r)xg^u + {r - q)i,g^x + {q - k)^g^x] + ig'^ig^s 
sm9wtg9w + g{oi2 - ai)sm.9w + 5'(«8 - a9)cos6'w - al3g?>(ic9w][q^,gux - qu 
gtix] + ig‘^{-gausm9wtg9w + ^gsec9w}[kf,gux - r^g^xl 
igsm9w[{k - r)xg^,y + (r - q)^,gyx + {q - k)ug^,x] + ig^{-a3Sin9w + (ai - 
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a2)cos6'vy + (as - ag)sm9w}[qf,gux - qug^^x] + ig^ausm9w[k^,gux 
-r^g^ix] (31) 

iii) Here, we write the Feynman rnles corresponding to the different vertices with 
4 gange electroweak bosons (Iyiy 2 V 3 V 4 )) shown in Fignre 


^w-w-w+w+ 

^Z0Z2W-W+ 

^A^Z 0 W^W+ 

^Ai,A,W^W+ 

^ZOZOZOZO 


+ fi'^(«4 - as) + 2 g‘^{a 3 + ag + ai 3 )[ 2 gf,^gxp - g^xgup - g^pgux] + 2ig'^ 

(a4 -|- ot3)\(]p\gi,p + g^pgux] 

ig‘^{-cos‘^9w - 2g‘^a3 + g‘^sec^9w{a^ + a^)}[2gpygxp - g^xgup - g^pgux] + *5'^ 
{(a4 + as + ae + a7)sec^6'ty + ai2tg^9w - ausm‘^9wtg^9w}[gpxgup + g^pgux] 
-ig^sm9wcos9w[l + g‘^sec^9wa3][2gpi,gxp - g^xg^p - g^pg^x] + ig*{sm^9w 

tg 6 *vyaii — -tg9w<yi2}[gpxgup + g^pgux] 

-ig‘^{sm‘^9w[2gp^gpX - g^xgup - g^pg^x] + g‘^si-^^Owaii[gpxgvp + giipgvx]} 
2isfic^9wg%OiA + as) + 2(a6 + ar + aio)] [g^vgxp + g^xgup + g^pgux] (32) 
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Table Captions 


Table In this Table we display the relations between different sets of chiral Lagrangian 
parameters. They correspond to those used by Longhitano |]T^ , Feruglio |T^ and Dobado 
et.ah p, in columns 1, 2 and 3, respectively. 


Table |^. In this table we show the sensitivity functions and associated to each 

ttfc {k = 0,1, 2, 3,4, 5) parameter. They were calculated with the minimal cuts and under 
the hypothesis of linear behaviour of c) with respect to each 


Table ^ Optimal cuts to calculate the statistical significance for typical positive 

and negative ak values {ak = ±0.005) in the or W^Z^ final states. The first entry 
in column 3 refers to the minimum invariant mass and the second to the minimum allowed 
Ptz, both them given in GeV. The total number of events in the Zero Model and with 
the typical chosen ak obtained after having applied the optimal cuts, are represented in 
columns 4 and 5. They correspond to 1 year of running for the LHC and a maximal 
rapidity of 2.5 for the final bosons. 

Table ^ Here, we shown a comparison between the statistical significance function 
given by the typical ak in Z^Z^ or W^Z^ final state, for the minimal and the optimal 
cuts. In both cases, we have fixed the maximal rapidity cut of 2.5 for two gauge final 
bosons, and 1 year of working time for the LHC. 


Table |^. In column 3, we display the statistical errors, (Aafc)stat., corresponding to each 
ak parameter in the Z^Z^ or channels. They were obtained with the optimal cuts 

(column 1) and under the assumption of linear dependence of the total number of events 
on ttfc. 


Table ^ Comparison of the total number of Z^Z^ or W^Z^ events obtained in the 
Zero Model applying the optimal cuts and choosing different structure functions: EHLQ, 
MRSD— and GRVHO. In the last three columns, we represent the statistical significance 
to estimate these imprecisions, 'rstruc.(c) eq.(|^), between EHLQ and MRSD— (column 
6 ), EHLQ and GRVHO (7 column) and MRSD— and GRVHO (8 column). These 
uncertainties are called, respectively, ri^., r^l and r^l- In all these computations we have 
taken |/imax = 2 / 2 max = 2.5 and an integrated luminosity for the LHC of 3 x 


Table |^. In this Table we represent the total number of or W^Z^ pairs obtained at 
LHC using the optimal cuts and setting j/imax = l/ 2 max = 2.5 and an integrated luminosity 
for the LHC of 3 X h In columns 2 and 3 we display the results in the Zero Model 
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and with constant g and g' conplings. Again in colnmns 4 and 5 we show the total nnmber 
of or W^Z^ events bnt in the Running Zero Model and with g{^/l) and 

eg. (p8[) ■ In the last two colnmns we display a measnre of the statistical signihcance that 
corresponds to inclnde the dependence on the energy of the parameters and conplings, 
given by the or fnnctions, eg. (p^) . 
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Figure Captions 


Figure Propagators of the tt^ and W^. 

Figure |^. Vertices with three gauge bosons or two gauge bosons and one GB, obtained 
from £nlsm eq.(0). 

Figure Different four gauge boson vertices that appear in the theory described by the 
lagrangian £nlsm eq.(0). 

Figure Different diagrams to one loop order contributing to gg —>■ cross section, 

in the SM without a Higgs boson. 

Figure We display the diagrams at tree level order that contribute to the cross section 
qq Z^Z^ in the framework of chiral perturbation theory (with Dnlsm)- In this case we 
obtain the same result as in the SM without a Higgs boson. 

Figure ^ We show the only diagram contributing to the different Z^Z^ Z^Z^ helicity 
amplitudes obtained with £nlsm eq. dD. at tree level order. 

Figure The scattering amplitudes for the process W^W~ —> calculated with 

-^NLSM eq.(^ at tree level order in the Landau gauge, receive contribution from all these 
channels. 

Figure ||. The cross sections qq' —> W^Z^ are calculated adding the contribution of the 
diagrams displayed in this Figure. They are obtained with £nlsm eq.(|^ at tree level 
order and in the Landau gauge. The initial quarks are supposed to be massless. 

Figure ^ Different diagrams contributing to the W^Z^ W^Z^ scattering amplitudes. 
They are calculated with Dnlsm eq.(0) at tree level order and in the Landau gauge. 

Figure |T^. The scattering amplitudes W^Z^ are obtained with the Lagrangian 

given in eq.(^) adding the contribution of the diagrams displayed in this Figure. The 
calculation is made at tree level order and hxing the Landau gauge. 

Figure |^. This graphic represents the sensitivity of the LHC to the ao parameter in 
and W^Z^ channels. We display the total number of Z^Z^ (solid line) and W^Z^ 
(dashed line) events obtained in the LHC, for different values of the ao parameter. The 
kinematical cuts are the minimal ones. The other ak parameters, apart from ao, have 
been set to zero. 
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Figure 12. This Figure represents the sensitivity function to ai parameter for the LHC, 
using the same kinematical cuts as those in Figure 0. We display the total (solid 

line) and (dashed line) events obtained versus ai values. 


Figure |T^ This Figure represents the sensitivity function to a 2 parameter for the LHC. 
As in Figures 0 and 0 the total number of and W^Z^ events versus a 2 values. 


Figure This Figure represents the sensitivity function to parameter for the LHC, 
using the same kinematical cuts as in previous Figures. We show the total (solid 

line) and (dashed line) events versus a^- 


Figure |^. This graphic represents the sensitivity of the LHC to the Q !4 parameter in the 
and W^Z^ channels. We display, like in Figures |^, |^, ^ and the total number 
of Z^Z^ (solid line) and (dashed line) events obtained at LHC, for different values 

of the 04 parameter, using the same kinematical cuts as in previous Figures. 


Figure [T^. Sensitivity function to parameter for the LHC. We have used the same 
kinematical cuts as in previous Figures. We display the total (solid line) and 

(dashed line) events obtained versus values, as we did in hgures 0 to pA[ 


Figure 17. In this Figure we study the variation of the total number of Z^Z^ (solid line) 
and W^Z^ (dashed line) events with the cut on the minimal invariant mass, We 

have set all the parameters to their values in the Zero Model and the cuts Przmin) 2/imax 
and |/2max to their minimal values. 


Figure As in Figure 0, we display the behaviour of the total number of Z^Z^ (solid 
line) and W^Z'^ (dashed line) events in the Zero Model versus the cut in the minimal 
transversal momentum of the Z^ boson. We have taken the other kinematical bounds as 
their minimal values. 


Figure |^. Here we reflect the dependence of the total number of Z^Z^ and W^Z^ events 
on the maximal rapidity bounds of hnal gauge bosons, |/imax = 2 / 2 max- We work in the 
Zero Model and £x the minimal cuts in invariant mass and transversal moment. The solid 
line corresponds to the Z^Z^ final state and the dashed one to the produced W^Z^ pairs. 


Figure ^ Statistical signihcance corresponding to the ao parameter. We display the 
ro function eq. (p!8[) for different values of ao with respect to the Zero Model. The solid 
line represents the Z^Z^ hnal state and the dashed line the W^Z^. We use the same 
kinematical cuts as in Figures 0 to 1^ 
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Figure Statistical significance corresponding to the ai parameter. We show as in 
Figure |^, the (solid) eq.(|^ and (dashed) functions with respect to the Zero 


Model versus these ai values. We use the same kinematical cuts as in Figures pd] to p! 6 . 


Figure This Figure represents the statistical signihcance for some values of a 2 with re¬ 
spect to the Zero Model. We display as in Figures 0 and 1 ^, and for the same kinematical 
cuts, the and eq.(|^ functions. 


Figure This Figure reflects the statistical signihcance of 03 . We show the eq.(P!B|) 
functions for Z^Z^ (solid) and W^Z^ (dashed) with respect to the Zero Model versus 
values, using the same kinematical cuts as in previous Figures. 


Figure |^. Statistical signihcance corresponding to the Q !4 parameter. As in Figures pO|to 
^ we show the (solid) (eq. |T^ and (dashed) with respect to the Zero Model 

versus the 04 values. We have used the same kinematical cuts as in previous hgures. 


Figure ^ The same as in Figures ^ to ^ for the parameter. 
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Table 1: 


tto 

O 

to 


CXl 

g'/gai 

LigtgOw 

Oi2 

g7ga2 

—LgR/2ig6w 

as 

—as 

—LgL/2 

C(4 

04 

L2 

«5 

0-5 

Li 

ae 

dQ 

— 

Oil 

dj 

— 

as 

—Os 

— 

ag 

—Og 

— 

aio 

2aio 

— 

ail 

Oil 

— 

«12 

2ai2 

— 

«13 

^13 

— 


Table 2: 


Oik 



ao 

-899.54 

2779.89 

ai 

1241.27 

-5408.20 

a2 

-1481.65 

6150.98 

as 

6751.24 

-48366.79 

0^4 

2445.43 

-18981.46 

as 

5924.90 

-20142.33 
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Table 3: 


ttk value 

Channel 

iiflfMffpiH 

iVh)(a°;c°P) 

iVh)(a;c°P) 

0:3 = —0.005 


( 200 , 10 ) 

16152.55 

16119.28 

as = —0.005 


(200,300) 

522.59 

581.41 

as = 0.005 

z^z^ 

( 200 , 10 ) 

16152.55 

16186.47 

as = 0.005 


( 200 , 200 ) 

2268.62 

2160.74 

a 4 = —0.005 

w^z^ 

(1150,500) 

32.29 

66.48 

a 4 = 0.005 

z^z^ 

(1150,400) 

13.50 

20.01 

as = —0.005 


(1150,400) 

61.56 

85.78 

as = 0.005 

z^z^ 

(1150,400) 

13.50 

28.70 


Table 4: 


Ok value 

Channel 

rk{oi-c^^) 

rfc(a;c°P) 

as = —0.005 

Z'^Z^ 

0.26 

0.26 

as = —0.005 

W^Z^> 

1.01 

2.57 

as = 0.005 

z'^z^ 

0.27 

0.27 

as = 0.005 

w^z^> 

0.51 

2.27 

a 4 = —0.005 

w^z^ 

0.45 

6.02 

a 4 = 0.005 

z'^z^ 

0.18 

1.77 

as = —0.005 

w^z^ 

0.39 

3.09 

as = 0.005 

z'^z^ 

0.37 

4.14 
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Table 5: 



Channel 

(^Clfc)stat. 

(200,300) 


(Aa3)stat. = 1.95 X 10-^ 

(200,200) 


(Aa3)stat. = 2.20 X 10“^ 

(1150,500) 

w ^ z ° 

(Aa4)stat. = 8.31 X 10-" 

(1150,400) 

z ^ z ^ 

(Aa4)stat. = 2.82 X 10“^ 

(1150,400) 


(Aa5)stat. = 1.62 X 10“^ 

(1150,400) 

z ^ z ^ 

(Aa5)3tat. = 1.21 X 10-^ 
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Table 6: 


Channel 

Optimal cuts 

iVh)(a°;c°P) 
EHLQ{set II) 

iV«(aO;c°P) 
MRS - D' 

iV(*)(a°;c°P) 

GRVHO 

'12 

.0) 

^13 

'23 

ZZ 

(200,10) 

16113.31 

17242.93 

16617.27 

8.90 

3.97 

4.76 

ITZ 

(200,200) 

2298.57 

2526.75 

2452.56 

4.76 

3.21 

1.48 

WZ 

(200,300) 

535.80 

594.34 

578.75 

2.53 

1.86 

0.64 

PTZ 

(1150,500) 

36.02 

40.79 

39.97 

0.79 

0.66 

0.13 

ZZ 

(1150,400) 

14.92 

15.88 

15.39 

0.25 

0.12 

0.12 

PTZ 

(1150,400) 

68.02 

76.72 

74.78 

1.05 

0.82 

0.22 


Table 7: 


Optimal cuts 

No running 

Running 

run. 

(a/ ^min) 






■nygn 

(200,10) 

16149.13 

91141.49 

15355.82 

86149.13 

6.24 

16.54 

(200,200) 

— 

2281.63 

— 

2093.83 

— 

3.93 

(200,300) 

— 

532.33 

— 

505.18 

— 

1.18 

(1150,500) 

— 

36.01 

— 

45.75 

— 

1.62 

(1150,400) 

14.91 

68.01 

11.28 

78.07 

0.94 

1.22 
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